On 
(N 



X 



DESY 07-198 
hep-ph/0711.2196 



Baryon scattering at high energies: 

C< ! wave function, impact factor, and gluon radiation 

O 

o 

>. 

O. J. Barters a) * and L. Motyka a) ' 6) 

a ) // Institute for Theoretical Physics, University of Hamburg, Germany 
SZ • b ^ Institute of Physics, Jagellonian University, Krakow, Poland 

g^ November 14, 2007 

43' 



> 

^O ! Abstract 



The scattering of a baryon consisting of three massive quarks is investigated in the high energy limit 
of perturbative QCD. A model of a relativistic proton- like wave function, dependent on valence quark 
longitudinal and transverse momenta and on quark helicities, is proposed, and we derive the baryon 
f ^ impact factors for two, three and four i-channel gluons. We find that the baryonic impact factor can be 

written as a sum of three pieces: in the first one a subsystem consisting of two of the three quarks behaves 
very much like the quark-antiquark pair in 7* scattering, whereas the third quark acts as a spectator. The 
second term belongs to the odderon, whereas in the third (C-even) piece all three quarks participate in the 
scattering. This term is new and has no analogue in 7* scattering. We also study the small x evolution of 
gluon radiation for each of these three terms. The first term follows the same pattern of gluon radiation 
as the 7*-initiated quark-antiquark dipole, and, in particular, it contains the BFKL evolution followed 
by the 2^4 transition vertex (triple Pomeron vertex). The odderon-term is described by the standard 
BKP evolution, and the baryon couples to both known odderon solutions, the Janik-Wosiek solution and 
the BLV solution. Finally, the t-channel evolution of the third term starts with a three reggeized gluon 
state which then, via a new 3-^4 transition vertex, couples to the four gluon (two-Pomeron) state. 
We briefly discuss a few consequences of these findings, in particular the pattern of unitarization of high 
energy baryon scattering amplitudes. 
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1 Introduction 

In recent years deep inelastic electron proton or electron nucleus scattering (DIS) at small x has attracted 
much interest, and it has stimulated intense studies of high energy QCD. At high energies, the total cross 
section of a virtual photon scattering on a target, in a first approximation, can be described in terms 
of a photon impact factor and a Balitsky-Fadin-Kuraev-Lipatov (BFKL) Green's function [1-3]. When 
restricting to the large N c limit, and assuming a large target, unitarity corrections to this first approximation 
are described by the nonlinear Balitsky-Kovchegov (BK) equation [4, 5] which, in the language of BFKL 
Green's functions, represents the infinite sum of fan diagrams [6]. The BK equation was initially obtained 
in the s-channel color dipole picture (in the large N c limit) [7,8]. Beyond the large N c limit one has to 
include the full color structure of the 2^4 reggeized gluon vertex [9] which leads to the Balitsky hierarchy 
of integral equations [4] or to the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) 
equations [10]. In many of these calculations the incoming virtual photon plays a vital role: its large 
virtuality Q 2 justifies the use of perturbation theory, and its impact factor consists of a quark-antiquark 
pair which forms a color dipole configuration. This simple structure is also intimately connected with the 
fan-like structure of the diagrams resumed by the nonlinear BK equation. 

The advent of the LHC challenges us with the task of developing a theoretical understanding of scattering 
in high energy proton-proton collisions, which is related to the structure of unitarity corrections in baryon- 
baryon scattering. In this paper we will perform a study of the high energy behavior of baryon scattering 
within perturbative QCD. It is clear that the problem of high energy nucleon scattering is much more 
complex than it was in the virtual photon case. First of all, in nucleon-nucleon or nucleon-nucleus scattering 
the incoming projectiles are nonperturbative, and the accuracy of perturbative calculations is not under 
good theoretical control. We shall circumvent this problem by studying a fictitious scattering process of a 
heavy and small baryonium system, in analogy to the heavy onium proposed as a test case for perturbative 
unitarity corrections in DIS [8]. For such processes the perturbative calculations provide reliable results. 
Next, the baryonium scattering is expected to differ significantly from the onium scattering. The main 
reason is the difference of the color structure: in contrast to the color dipole the baryon is a color singlet 
formed by three valence quarks. Also, the application of the large iV c limit which played the crucial role 
in the construction of the dipole model is rather difficult in the baryon case: one needs exactly N c quarks 
to build the color singlet of the SU(N C ) group, and this system becomes rather complex for N c — ► oo. In 
fact, a few years ago, it was explicitly pointed out [11] that the simple picture of gluon radiation which 
has emerged in the QCD dipole picture does not work in the case of an incoming three quark color singlet 
system; however, no alternative solution had been derived. Thus, we shall address the issue of gluon radiation 
from three quarks at N c = 3, within a perturbative baryonic system and compare with the perturbative 
quark-antiquark system. 

The basic and universal object that characterizes properties of the baryon is its wave function. Inspired 
by the success of the concept of the photon wave function [7] which turned out to be very fruitful in studies of 
high energy scattering, we start from a local three-fermion quark current operator with the quantum numbers 
of the proton and construct a relativistic invariant infinite momentum frame wave function for the lowest 



Fock component of the baryon, consisting of three valence quarks. The resulting wave function contains 
a non-trivial dependence on quark helicities and angular momenta. For the current operator we chose the 
baryonic operator proposed by Ioffe [12], which has been shown to provide a reasonable phenomenological 
prescription of the nucleon properties [13]. In order to take into account the nonperturbative nature of the 
baryon we make use of the Borel transform technique which has been developed in the context of QCD sum 
rules. 

This paper is not intended yet to deal with a detailed phenomenology of the baryon structure and 
scattering, - - thus we do not attempt, for example, to tune the obtained wave function to describe the 
existing data on proton form- factors and high energy scattering. Nevertheless, apart from developing a 
theoretical laboratory for studying scattering of baryon states at high energies, one may hope that our 
perturbative analysis finds structures which remain also relevant beyond the perturbatively safe region. An 
extrapolation of our results on the heavy baryonium to the realistic proton case may, therefore, very well 
allow for some useful phenomenology. More detailed studies in this direction will be left for future work. 

Starting from integrals over squares of these baryonic wave functions and coupling t-channel gluons to the 
three quark lines we define baryonic impact factors, in close analogy with the photon impact factor in deep 
inelastic electron proton scattering. The small- x evolution of baryon scattering amplitude will be analyzed, 
again, following the strategy developed in the context of the virtual photon scattering [9, 14, 15]. First we 
consider, in lowest order, the elastic scattering of the baryonic system on a single quark: by coupling two 
t-channel gluons to the three-quark system, the baryonic impact factor is obtained. Three or four t-channel 
gluons appear if one considers, again at lowest order, multi-particle amplitudes, e.g. 3^3 processes in 
a suitably defined high energy limit. In the next step, one considers higher order diagrams in the leading 
logarithmic approximation: this leads to rapidity evolution equations, describing the radiation of gluons 
from the three-quark system. 

Our main results are the following. We propose a model of the baryon wave function with a non-trivial 
quark helicity and angular momentum structure. Then we express the baryon impact factor in terms of the 
wave function, for an arbitrary number of coupling gluons. The obtained baryonic impact factor can be 
written as a sum of several pieces, each of them having its own evolution equation. First, there is a term in 
which one pair out of the three quarks scatters whereas the third quark acts as a spectator. Although the 
two quarks which participate in the interaction are in a color anti-triplet configuration, they behave very 
much like the quark-antiquark pair in the photon case. In the lowest order, two t-channel gluons couple to 
this quark pair. In higher order the two gluons start to reggeize and to produce the full BFKL ladder, while 
the third quark of the baryon state remains an inactive spectator. Also, the well-known 2^4 gluon vertex 
appears, indicating the beginning of the same fan-like structure as in the quark-antiquark case. Altogether, 
this piece of the baryon impact factor radiates gluons in very much the same way as the quark-antiquark 
pair in the photon case. 

Next, there is the odderon term, similar to the one discussed in [16]: here all three quarks participate, 
and the t channel state carries C = — . In lowest order, three gluons couple to the three quarks; in higher 
order the state evolves according to the Bartels-Kwiecihski-Praszalowicz (BKP) evolution equation [17, 18]. 
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Figure 1: Energy discontinuity of the 2^2 process: 7* + q — ► 7* + q. 

Finally, a third, C-even, piece of the baryonic impact factor appears in which again all three quarks 
participate. This piece has no counterpart in the quark-antiquark case and, together with the odderon, it 
makes the baryon really behaving differently from the photon (or the vector meson). The state consists of 
one reggeized gluon with even signature and two usual odd reggeized gluons. It obeys the BKP evolution 
in the three Reggeon channel and it decays into four reggeized gluons via a new gauge invariant 3^4 
reggeized gluon vertex. 

The paper is organized as follows. We begin with a short section describing the general framework 
in which our calculations are carried out. We then (Section 3) turn to the baryon wave function which 
enters the baryon impact factor. In the following Section 4 we describe the baryon impact factor and its 
decomposition into the three pieces described above, and in Section 5 we discuss the rapidity evolution of 
these pieces. Section 6 contains a short discussion of the baryonic impact factor in configuration space, 
and in Section 7 we analyze the quark-diquark limit of the baryon wave function. Finally, in Section 8 we 
summarize our results and discuss a few potential implications. 

2 The framework 



In our calculation we will follow the analysis of the scattering of a virtual photon described in [14, 15]. In 
leading order the scattering of a virtual photon off a quark is described by the exchange of two gluons. The 
coupling to the photon is described by the photon impact factor, Z?2;0) which most easily is obtained by 
the energy discontinuity of a closed quark loop (Fig. 1). Making use of the Regge factorization, the same 
impact factor can also be used in other elastic scattering processes, e.g. in the scattering of a virtual photon 
on a heavy onium target. Higher order corrections, in the leading logarithmic approximation, lead to the 
reggeization of the t-channel gluons and to the exchange of a BFKL Pomeron between the photon impact 
factor and the target. 

If one is looking for corrections containing more than two reggeized t-channel gluons one has to go beyond 
the leading logarithmic approximation. In the elastic scattering process 7* + q — ► 7* + q , both leading 





Figure 2: Multiple energy discontinuities of the 3^3 process: 7* + q + q 
(a) lowest order diagram, (b) two examples of higher order diagrams. 
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order and NLO corrections retain the structure of a single ladder. A i-channel state with four reggeized 
gluons appears first in NNLO. A convenient way to avoid the complications connected with such a high 
order calculation is the study of multi-particle processes, e.g the 3^3 process 7* + q + q ^ 7* + q + q , 
the scattering of a virtual photon on two independent quarks (Fig. 2) in the triple Regge limit. This 
process depends upon three independent energy variables, and the triple energy discontinuity can be easily 
computed in the approximation where, in each order perturbation theory, the maximal number of large 
energy logarithms is kept. The lowest order contribution is described by the exchange of four gluons. In 
higher order, these i-channel gluons reggeize and start to interact. As discussed in detail in [14, 15], the 
all-order result can be cast into the two sets of diagrams shown in Fig. 3. 

The first term starts, at the photon impact factor, with a BFKL Green's function, then undergoes the 
transition into the four gluons and continues with the BKP evolution of the four gluon state. In the large- N c 
limit, the four gluon state turns into two noninteracting BFKL systems, i.e. we see the beginning of the 
fan-diagram structure of the BK equation. The second term consists of a simple BFKL Green's function, 
with higher order splittings of the reggeized gluons at the lower end. As a remarkable feature of this results, 
in both contributions only two reggeized gluons couple to the photon impact factor, despite the fact that 
diagrams with four gluons — such as the one shown in Fig. 2a — are included: the apparent 'disappearance' 
of these contributions is a result of the gluon reggeization which manifests itself in generalized bootstrap 
relations. 

The same strategy can be used to investigate t-channel states with higher number of t-channel reggeized 
gluons. For example, six gluons appear in the 8-point amplitude j* + q + q + q^j* + q + q + q, i.e. the 
scattering of a virtual photon on three independent quarks. The analysis of this case has been investigated 
in [15]. 

Although these results are — initially — derived in the context of a higher order multi-particle processes 
(e.g. the 3^3 scattering process), they nevertheless can be used also in a 2 — > 2 process. The diagrams 





Figure 3: Decomposition of the sum of all diagrams in Fig. 2b into (a) irreducible and (b) reggeizing pieces. 




Figure 4: Four gluon contribution to the Reggeon unitarity equation of elastic 7*7* scattering. 



shown in Fig. 3 satisfy Reggeon unitarity equations in all three t-channels. Taking the discontinuity across 
the four Reggeon state, the partial wave above this can be used to construct the four Reggeon state in the 
2^2 process shown in Fig. 4. 

In this paper we will apply the same construction, replacing the virtual photon by a three quark system. 
Modeled by the four fermion operator introduced by Ioffe in the context of the QCD sum rules [12, 19], 
the incoming 'baryon' splits into three quarks which then couple to 2, 3, or 4 gluons. In order to take into 
account the non-local nature of the incoming baryonic bound state we introduce a form factor: we employ a 
technique used in the QCD sum rules [19] and use the Borel transform of the perturbative expression [20, 21]. 
The exponential nature of this form factor also guarantees the convergence of the momentum integrals inside 
the impact factor. 




Figure 5: Multiple discontinuity of the impact factor for elastic baryon scattering. 

3 The baryon wave function 



We consider the multiple discontinuity of a non-forward baryon impact factor in elastic high energy scat- 
tering. Large momenta are directed along the z-axis, and the incoming and outgoing baryons move at 
small angles with respect to the z-axis, as shown in Fig. 5. Their momenta P, P' have a large "+" light- 
cone component, P + , and their transverse momenta are denoted by P, P', respectively. We introduce the 
light-like vector q^ = (q, 0,0, — q), q 2 = with s = (P + q) 2 ~ 2P • q, and we assume that s is large: 
s S> M 2 , P 2 , P' . The quark momenta pi are 
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For the longitudinal quark momenta it will sometimes be convenient to use the notation 



Pi 



a,P^ 



Pi 



A 



(2) 



We shall use p = 7 M p M = 7 • p for contraction of four- vectors and Dirac 7 matrices. The adopted model of 
the proton state is defined by 

(0|r?(0)|iV(P,A)) =A N w x (P), (3) 



where w\(P) is the proton spinor with momentum P and helicity A, 

v(x) = e K1K2K3 [(u^(x)) T C^u K2 (x)] W d K *(x) 



(4) 



is the baryonic Ioffe current [12], C is the charge conjugation matrix, and Ki are color indices. The Ioffc 
operator is not the only possible choice of the baryon current, — in the context of distribution amplitudes, 
the possible baryonic operators for the proton were classified in Ref. [22], and it was shown that the Ioffe 
current gives a rather good description of baryon form- factors [13]. We therefore chose, as a test case, the 
Ioffe operator to model the baryonic impact factor. 1 



x It is worthwhile to stress that our baryon wave functions are different from the distribution amplitudes. In the collinear 
approach one probes the baryon with a hard external scale, Q 2 , and the baryon structure is represented by series of distribution 
amplitudes with increasing twist, that is with increasing power- like suppression at large Q 2 . The distribution amplitudes depend 
on the quark longitudinal momenta, and they obey evolution equations in logQ 2 . In contrast to that, we are interested in the 
baryon wave function with full momentum dependence probed at a moderate momentum scale, and the evolution applies to the 
rapidity of gluons radiated from the baryon impact factor. 
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Figure 6: The proton vertex as given by the Ioffe current. 



For the calculation of the baryonic impact factor we will need the matrix elements (Fig. 6) in the helicity 
basis, 

[dxMl^wx(P)] ■ [^(pOVCVi^fo)]. (5) 



In the second term we can also write: 

[«Ai(Pl)y VA 2 (P2)] 

where v (in the Dirac notation) denotes the -u-spinor of the u quark. 



(6) 



3.1 The massless quark case 

Using the calculus described by Brodsky and Lepage [23] we compute the Dirac spinor matrix elements. 
The details of the calculations are described in Appendix A. For simplicity, we start from the massless quark 
case, and the case of massive quarks will be analyzed afterwards. Thus we obtain: 



[d\(p3)l5j»w x (P)] ■ [u Xl (pi) ^ Cj° u* X2 (p 2 ) 
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(8) 



2M <5_ Al Aa <^ 8 Xl ,xV\ ■ + ^A 2 ,A»?A • 



where the transverse complex vector rf X is defined by 



r tx = (l,i\), A = ±l, 



(9) 



and P = P1+P2+P3 is the transverse momentum of the incoming baryon. The cross product of two transverse 
vectors p\ = (pi,Pi) and P2 = (P21P2) should be understood as a number Pi x P2 = p\p\ — p\p%- It 
turns out that formula (7) may be re-expressed in a more compact form, by using the vectors r}\ with the 
following identity: 

(Pi-Vx) ipr-n-x) = (K + ^)(p|-iAp|) = ptpl+plvl + iKplpl-pipl) = PvP2 - i \P1XP2, (10) 
which holds for any pair of transverse vectors, p\ and P2- Using this relation one gets: 

[d\(P3)l5l»W X {P)] ■ [«Ai(Pl) 7 M C7 «Aa(P2)] 
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v /aia 2 «3 



«2 






n a 
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(11) 



In what follows, we shall express all formulae in this compact notation. 

Next we couple a gluon of momentum k = (3q + A; to one of the quark lines with momentum pi (Fig. 7) . 
Fixing the momenta of the outgoing quarks at pi, P2, and p%, the quark line to the left of the gluon vertex 
carries momentum pi — k. Using, at the gluon vertex, the eikonal approximation, one arrives at the spinorial 
factor q. With 2pi ■ q = ais 3> pf,k 2 , etc. , one obtains, for the upper u quark, 



u(Pi) Q (Pi ~k) = 2p t -q u{pi -k) + ... , 



(12) 



where . . . stands for terms which are power suppressed in s. An analogous expression holds for the d quark, 
whereas for the second u quark we use: 



(P2 ~ k) qu*(p 2 ) = 2p 2 -qu*(p2 - k) + . . . . 



(13) 



As a result, on the r.h.s. of Eqs. (12) and (13), the transverse momentum of the quark spinor coincides with 
the transverse momentum of the internal quark line next to the baryon vertex. The sum of the outgoing 
transverse momenta equals 

Pi+P2+Ps=P + k. (14) 

Matrix elements corresponding to multi-gluon couplings to spinor lines may be simplified by iterating 
Eq. (12) in the following way: 

u(p) q [7 • (p - ki)j q ... q[j ■ (p-ki- ...- k n )\ ~ (2p ■ q) n u(p - ki - . . . - k n ). (15) 

For completeness, we remind that, in the case of an outgoing antiquark, an additional minus sign appears: 

— (p — k) qv(p) = —2p-qv{p—k) + . . . , (16) 

This minus sign is due to the opposite direction of the momentum along the antifermion line. Similarly: 

[-7 • (p ~ fa - ■ ■ ■ - kn) ] q . . . g[-7 • (p - fa)} q v(p) ~ {-2p ■ q) n v(p - fa k n ). (17) 
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Figure 7: Single gluon coupling to the baryon — ► X transition. 



This change in sign plays a crucial role in the photon impact factor [7]. 

We are now ready to describe the amplitude for the process: baryon + gluon — ► 3 quarks, defined by 
the diagrams shown in Fig. 7. We define the shifted momentum of the upper quark 



pi =Pi-k, 



with 



and use 




(18) 
(19) 



CKl Oil OLJ, I 

(and analogous expressions for the gluon coupling to quark lines 2 and 3). We introduce the amplitudes 
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(21) 
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«2 "3/ " \«1 «3 



(22) 



Here the upper three indices of denote the helicities of the upper two u quarks with momenta p\ and 
P2, and the lower d quark with momentum P3, respectively. The subscript refers to the helicity A of the 
incoming baryon. We leave the normalization constant Mq unspecified here; the normalization will be 
fixed at the level of baryon wave function. The amplitudes for the diagrams shown in Figs. 7a- 7c are then 
simply obtained from (21), (22) by the replacements p\ — ► p[, P2, — ► p' 2 and P3 — ► p 3 , respectively. 2 



2 To be precise, the functions O give the momentum dependent part of the scattering amplitudes, up to a global normalization 
factor, that is proportional to the strong coupling constant g. Obviously, the color factors are not accounted for in (21) and 
(22), — they will be treated explicitly later on. 
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Note that, for each of the three diagrams, the denominator is just the energy denominator in non- 
covariant perturbation theory, for instance one obtains for Fig. 7a, 



Eh. 



aryon 



E 



3 quark 



±U 2 + P 2 




(23) 



The amplitudes should be invariant under Lorentz boosts in the transverse directions, parametrized by a 
four-velocity u^ ~ (l,it,0), \u\ <C 1: 



Pi 



Pi 



1 , 
Pi + -P>, 



Pi 



Pi - 



Pi 



(24) 



The numerators are manifestly invariant under these transformations, and the denominators may be also 
rewritten in an explicitly invariant form using the identity 
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p2 , Pi , P2 I ^3 
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(p 2 - a 2 P) 2 , (p 3 - a 3 P) 
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(25) 



(26) 



CKi «2 CK3 CKi CK2 «3 

The denominators have poles for the invariant mass of the three-quark system equal to the proton 
transverse mass. Clearly, this is a consequence of using a point-like vertex for the proton-quark coupling 
and neglecting the bound state effects. These effects cannot be described within perturbative QCD and 
should be modeled. Following [19, 20] we propose a model that preserves Lorentz and helicity structures of 
the perturbative expressions, where the bound state effects are absorbed into the Borel transform. 

The Borel transform of a function f(s) is defined in the standard way: 

■n+l / d N" 



B s [f](M 2 B ) = lim 



n^oo n\ 



*> m - 



CO, 



s/n -> M 



B; 



(27) 



where Mb is the Borel parameter. In order to model the baryon scattering amplitude we shall apply two 
independent Borel transforms w.r.t. the negative virtualities: Q 2 = —P 2 of the incoming and Q' 2 = —P' 2 
of the outgoing baryon, to the perturbative amplitudes obtained with the point-like vertex. Formulae (21) 
and (22) were presented for P 2 = M 2 . The corresponding formulae for general virtualities are obtained by 
substitutions M 2 — ► P 2 in the denominators. In the baryon impact factor, the virtuality P 2 appears only in 
the energy denominator of the vertex amplitude K 1 ' 3 ({aj}, {pi}; P) of the incoming baryon, and the 



*(Ai,A 2 )A 3 , 
A 



er" a " V3 ({ai},te};P / ; 



virtuality P' 2 only in the denominator of the amplitude 

(see Sec. 4 for more details). Therefore the two Borel transforms may be performec 

0, that is already at the level of the baryon wave function. Thus we evaluate 

1 



B f 



Q 2 + M 2 X 



(M 



exp (-M 2 X /M 2 B ) . 



of the outgoing state 
independently for each 



(28) 



This result, applied to the amplitudes 0, leads to the substitution 



exp 



P 2 +P 2 -J2 



3 

i=l a 



V\ 




(29) 



11 



Before we complete the model we shall perform some simplifications. We shall absorb into the wave 
functions a phase space factor (aia^as) -1 that appears in the baryon impact factor as a result of on-mass- 
shell conditions of the cut quark lines. In this way, the factor y/ a.ia.2013 present in the amplitudes Q will be 
removed from the wave functions. Obviously, the integration measure will be suitably modified as well. For 
simplicity, we introduce a normalization constant, M, of the wave function that will be fixed later. Thus, we 
choose the natural value of the Borel parameter Mb = M and obtain a model of the baryon wave function, 
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(30) 



x M \6 XuX r, x .(»-*) + S^to-l*-*) . (31) 

Clearly, the functions ^ given by Eqs. (30) and (31) are symmetric under the interchange of the u quarks, 
labeled by 1 and 2. When combined with the anti-symmetry in the color degrees of freedom it implies 
that the full wave function is anti-symmetric under interchange of the u quarks, as it must be. Interestingly 
enough, a similar Gaussian form of the wave function was proposed long ago [24] and it was shown to provide 
a good description of the nucleon form-factor data [25, 26]. An important difference of our model, however, 
is the presence of angular momenta of the quarks. The baryon angular momentum structure following from 
the model is most transparent in the coordinate representation and will be discussed is Section 6. 

The above derivation of the baryon wave function is based on perturbative QCD methods combined 
with the Borel transform technique. Clearly, we are not able to control the accuracy of this procedure for 
the proton as it is a genuine non-perturbative object. Therefore the obtained wave functions can be only 
considered as a theoretically inspired model of the proton wave function. Therefore, in the next part we give 
the formulae for the wave function of a baryon consisting of three quarks with the same mass m, coming 
in two different flavors. These formulae will permit to consider the fictitious case of a large quark mass, 
for which the baryon becomes heavy and small, and the perturbative computation of its wave function and 
scattering is formally justified. 

3.2 Massive quarks 

We now apply the procedure described in the previous section to the case of the massive quarks. We skip 
the details of the derivation and present the result for the helicity amplitudes G of the transition: baryon 
to quarks, in which all three quarks were assumed to have the mass to: 
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(32) 



Note, that using relation (10) and taking m — ► one easily recovers formulae (21) and (22). The above 
formulae are promoted to the baryon wave functions *!', by going through the same steps as in the massless 
case, i.e. using the Borel transform and absorbing the a factors into the phase space factor. The final 
expressions for the wave functions are obtained from Eq. (32) by the replacement which combines both 
steps: 
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(33) 



4 Baryon impact factors 

4.1 General structure 

The amplitudes ^ may be combined with their complex conjugates to obtain the baryon impact factor. 
For the case of two gluons coupled to lines 3 and 1 we illustrate one example in Fig. 8. It was shown in 
the previous section that, in the high energy limit the spinorial part of the multiple discontinuity can be 
expressed in terms of universal matrix elements given by Eqs. (7) and (8), where the momenta of the quarks 
are evaluated at the quark-proton vertex. Also the denominator is determined by the virtuality of the quark 
to which the first gluon couples, and the virtuality can expressed in terms of the momenta of the quarks at 
the proton vertex. Thus, the impact factor can be obtained from overlap integrals, i.e. products of wave 
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Figure 8: A contribution to the two gluon exchange in baryon-quark scattering. 

functions with suitably adjusted momenta. As an example, we specify the overlap function corresponding 
Fig. 8: 



T xx '({k,k'};P,P') 



Ai,A2,A3 



[d 2 pi] [dcti] 



$ (Ai,A 2 )A 3 

A 



({aJ,{pi,p 2 ,p 3 };F') 



X ^ 



(Ai,A 2 )A 3 
A 



{{(Xi},{pi,P2,P 3 };P), 



with the integration measure 

[d 2 p t ] 



d 2 p\ d 2 p2 d 2 p3 



(34) 



(35) 



[dcci] = da\ dct2 da% 5{a\ + a 2 + «3 — 1). 

In the overlap functions J 7 , the upper helicity labels refer to the incoming and outgoing baryon states, 
respectively. Analogous overlap functions are defined for the other gluon couplings, and for the full impact 
factor we will have to sum over all diagrams. When evaluating the sum over the intermediate helicities 
Ai, A 2 , and A3 and summing over all diagrams, one finds, for the forward direction P = P' = 0, helicity 
conservation, i.e. the impact factor vanishes for A = —A'. 

Before including the remaining energy integrals and the color factors we generalize to the case of 3 and 4 
t-channel gluons. As outlined in Section 2, we have to consider multiple energy discontinuities. An example 
is shown in Fig. 9. For one of the discontinuity lines (in the case of Fig. 9, the central line) we fix the 
intermediate quark momenta and denote them bypi, p 2 >P3- The corresponding overlap function (Fig. 10) 
takes the form 



t AV ({ki, k'^PiP') = Yl 

Ai,A2,A3 



\d 2 pi] [doii 



^(Ai,A 2 )A 3 

A 



{{ai}, {pi + fc;,P2 + fc' 2 ,P3 + k' 3 };P') 



,(Ai,A 2 )A 3 
A 



x ^1.^-v, ({ ai }, { Pl - *i,p 2 - fc 2 ,p 3 - k 3 };P) , (36) 



where ki (k'A is the sum of momenta delivered by the gluons to the quark line i to the left (right) of the 
central cutting line. Obviously, 

J> + J>;.=F'-F. (37) 
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Figure 9: A contribution to the four gluon exchange in the scattering of a baryon on two independent quarks. 



p. -k. 




Figure 10: Example of a diagram defining the overlap function for a multiple discontinuity in Fig. 9. 
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Figure 11: Baryon impact factor. 

Since, as a consequence of the exponential form factor, all transverse momentum integrals are finite, we 

are allowed to shift, for each diagram separately, the loop momenta, such that to the right of the incoming 

baryon vertex, the momenta become pi, P2, andj>3. At the outgoing baryon vertex the momenta arepi + Zi, 

P2 + h, and p^ + h, where U is the sum of transverse momenta of all gluons coupled to the quark line i 

(Fig. 10): 

U = Y.ikj + Vj). (38) 

jeLi 



In the general case the overlap function can be written as: 



F xx '({h};P,P>) = Y, j\d 2 Pi][dai] [* ( i 1 ' Aa)A3 ({ai},{>i + Zi};i ,, )J v& ( A Al ' A2)Aa ({a t }, {p t }; P) . 

Ai,A2,A3 

(39) 

We now complete the definition of the baryonic impact factor (see Fig. 11). The general impact factor 
1\b for the transition A —> B , with TV gluons being coupled in the eikonal approximation, is defined as 



(_Y) _ J djhdfh djhj-i g££i q^ 2 ■ • • g m - n; „ c/ij / !j ../ x 

2tt 2tt '" 2vr s 



'AB 



DiscAr_i S^ 



(40) 



where S^ 12 '"^ N represents the amputated transition amplitude. In particular, for the elastic scattering of 
a single quark inside the baryon one obtains: 



I? = ^ Tr[t a »t a »-i...t a i] /(JO, 



N, 



where 



4f = (-ig) N ■ 



The baryon impact factor for TV t-channel gluons is then given by 

B xx '({l t };P,P') = /W Y, F XX '({h};P,P') C w (diagram) 



(41) 
(42) 

(43) 



diagrams 



where the color factor reads: 



Cn (diagram) 



3! 



[w- 1 ...* 01 ]^ 



j.bm-1-bm—i A)\ 



[fn^- 1 ...t cl ] K , 



3«3 



(44) 



In (43) the sum extends over all diagrams, and the numbers I, m and n of gluons, that couple to quark 
lines 1, 2 and 3, take all possible values between 1 and iV with the constraint I + m + n = N. The overlap 
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a) b) 

Figure 12: Diagrams defining B2fl. 

function J- xx ({li};P,P') is evaluated for each diagram separately, and in each case it contains a global 
delta function of the transverse momenta: 



F xx '({i,i'y,p,p') = F xx '(i 1 ,i 2 ,i 3 ) 5^ fj^ + p-pM 



We impose the normalization condition: 

F XX '(hJ 2 ,l3) 



rAA' 



h=h=h=o 
Correspondingly, we also extract a delta function from the impact factor: 



^'({Z,};P,P') = B xx \{l i })6w(j2li+P-P'\ 



(45) 



(46) 



(47) 



In the following we will restrict ourselves to the forward direction, P = P' = 0. Because of helicity 
conservation for the impact factor, we always have A = A', and we will drop the upper helicity labels, i.e. 
F aa '(ZiJ 2 J3) -► F(h,l2,h) etc. We will go through the cases ot N = 2, N = 3, and N = 4 gluons. We 
therefore define, for fixed N, the functions B^-q ({h}) projected on the C-even channel through 

B XX ({U}-P,P') \° p ~ = ™2 = B N ., ({li}) 6™ feli) , (48) 

and analogously for the C-odd projections i?/v;0 • 

In the remaining part of this section the main emphasis will be on the color structure of the impact 
factors, contained in Eq. (44). As the main result, we will find a decomposition into a sum of terms which, 
as it will be demonstrated in the subsequent section, stays invariant under evolution in rapidity. We stress 
that the results which follow are valid for an arbitrary overlap function F(li,l2,ls) , i.e. they do not rely 
on a particular model of the baryon wave function, provided the baryon has the valence degrees of freedom 
of three quarks. 



4.2 Two gluons 

We begin with the two-gluon coupling which is C-even. All diagrams are proportional the color tensor <5 aia2 , 
and it is suggestive to group them into three sets: in the first one, the two gluons couple to the quark pair 
(12), and quark 3 acts as a spectator (Fig. 12). In the second one the gluons couple to (13) and quark 2 
acts as spectator etc. Inside each set, we have four terms. We thus find: 



jU.2}/ 



,{1.3}, 



,{2,3}, 



5 2 ;0(fcl,*2) = 5 aia2 D^(k 1 ,k2) + Dl^ t (k 1 ,k2) + D^ t (k 1 ,k 2 ) 



(49) 
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a) b) c) 

Figure 13: Examples of diagrams defining B^-q and B^-q. 



with 



£> 2 { Jj 2} (fci,fc2) 
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2;0 

{2,3} 
2:0 



(fcl,fc 2 
(k 2 ,k 3 



12 



12 
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[F(k, 0, 0) + F(0, k, 0) - F(fci, fc 2 , 0) - F(fc 2 , fci, 0)], 
[F(k, 0, 0) + F(0, 0, fc) - F(fci, 0, fc 2 ) - F(fc 2 , 0, fci)], 
[F(0,k,0) + F(0,0,k)-F(0,k 1 ,k 2 )-F(0,k 2 ,k 1 )], 



(50) 

(51) 
(52) 



where fci, k 2 denote the gluon momenta and k = k± + k 2 . On the r.h.s. of (50) — (52), the momentum 
arguments of the F functions indicate which diagrams they represent: in the first (second) term of (50), 
both gluons couple to quark line 1 (2). In the third term, the first gluon couples to line 1, the second 
to line 2, and so on. The relative signs arise from the color structure. As a striking result, on the r.h.s. 
of (50) — (52), in each line the four terms have the same structure as the impact factor of the photon. In 
particular, each set satisfies the Ward identities, i.e. it vanishes as any of its momenta goes to zero. Since 
the pair of scattering quarks {i,j} is in a color anti-triplet state, one might, at first sight, interpret this 
set as the elastic scattering of an 'anti-triplet dipole'. However, it is important to stress that these three 
dipole-like components -D 2 . , are not independent from each other: the diagrams where two gluons couple 
to the same quark line, say, line 3 in Fig. 12a, contribute both to the pair (13) and (23). In this sense, 
one better views these quark pairs as 'anti-triplets inside the baryon'. Also, these configurations where one 
quark pair interacts whereas the third quark remains a spectator, should not simply be viewed as 'diquark 
states': in transverse coordinate space, the spectator quark can be far away from the quark pair (see the 
discussion in Section 7). One should also add that the normalization of the dipole-like components -D^o of 
the baryon impact factor is exactly 1/2 of the normalization of the genuine color dipole impact factor. At 
the two-gluon level, our results coincide with results of Ref. [27]. 

If, instead of our model for the baryonic impact factor, we would have used a completely symmetric 
baryon form-factor F^ s ' (which does not discriminate between u and d quarks) we would have arrived at a 
familiar result [28]: 



,(«) 



£$(*!, *2 



-9 



S ai 



FW(fc,0,0)-FW(fc 1 ,fc 2 ,0) 



(53) 



4.3 Three gluons 

In the case of three gluons (Fig. 13) we have to distinguish between even and odd C parity: in the color 
trace Eq. (44) we find both color tensors, f a ^ a ^ a i and d aia2a ' a . The first one belongs to even (Pomeron), the 
second to odd (odderon) C parity. 
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The C-odd baryon impact factor reads [16, 29, 30] 

S 3 ;0(*l,*2,* 3 ) = d™" 3 E 3 .n(k 1 , k 2 ,k 3 ), 



(54) 



where 



F 3 -o(ki,k 2 ,k 3 



24 



E 



(55) 



(56) 



2F a {k 1 ,k 2 ,k 3 ) -^F a (fci,fc -fci,0) + F a {k, 0,0) 

i=l 
and -F CT denotes the i* 1 functions, with its arguments being permuted by the permutation a: 

F <T (h,h,h) = F (I a (1) J a (2) J a (3)) ■ 

In E 3 -q the t-channel three gluon state is Bose symmetric, that is 

^3;o(*«r(l)»*<T(2).fc<r(3)) = E 3-,o(kl, *2, *3) (57) 

for any permutation a, and it obeys the Ward identities: 

E 3 -o(ki,k 2 ,k 3 ) = for any fcj -► 0. (58) 

On the r.h.s. of Eq. (55), the momentum structure of first term indicates that the three gluons couple to 
three quarks. The second and third term play the role of subtractions. This leads to the interpretation 
that, in this piece of the baryonic impact factor, in contrast to the structure found previously for 2 gluons, 
all three quarks participate in the interaction. Since each of the three gluon has negative C parity, this 
t-channel belongs to the C = — (odderon) state. 

For a completely symmetric model for the baryon form- factor expression (55), again, reduces to a known 
result [29]: 

3 

2F«(* 1 ,fc 2 ,fc 3 ) -Yl Fi3) (ki,k - k i: 0) + F^(k, 0,0) 



B$(fcl,* 2 ,fc 3 ) = ^d<™ 3 



(59) 



Next we turn to the terms proportional to f aia2Cl ' A which turn out to belong to even C . They can be 
grouped in the same 'dipole-like' form as in (49): 



B3-o{ki,k 2 ,kz) = D^ 2} (k 1 ,k 2 ,k 3 ) + D^,f } (ki,k 2 , k 3 ) + D^ 3} (k 1 ,k 2 ,k 3 ), 



where the dipole-like components have the structure known from the photon case, 



D^faMM) 



_ „ tai^a-i 



>{ij}, 



\{hj} 



{iJ}, 



Dffiik! + k 2 , k 3 ) - D$* (An + fc 3 , k 2 ) + D™'(k 2 + fc 3 , ki 



(60) 



(61) 



As in the photon case, the argument structure indicates the beginning of the reggeization of the gluons: for 
example, in the first term, the first two gluons with momenta k\ and k 2 'collapse' into a single reggeized 
gluon with momentum k± +k 2 . The t-channel system thus consists of two reggeized gluons only and hence 
belongs to C-even. In the next section we will show that this structure is preserved in the rapidity evolution. 

In the following it will be convenient to use a shorthand notation by writing, instead of D^q (ki+k 2 ,k 3 ), 
simply Z)£,f } (12, 3) etc. 
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4.4 Four gluons 

In the case of four gluons (Fig. 14) the color trace (44) contains //, dd, fd, and 55 color tensor structures. 
Beginning with the fd pieces, we find that they can be expressed in terms of the E'-function (55) which we 
have obtained for the odderon channel: 



£ 4; o(l,2,3,4) 



(62) 



- (f aia2b d ba3ai E 3 . (12,3,4) + / aia36 d ba2a4 .E3 ; o(13,2,4) + f aiaib d ba2a:i E 3 . (U, 2, 3) 
+ f a2asb d baia4 E 3 . (23, 1, 4) + f a2aib d baia3 E 3 . (24, 1, 3) + f a3Cl4b d baia ' 2 £ 3;0 (34, 1, 2) 



We then interpret this contribution as the odderon configuration with one reggeizing gluon. It agrees with 
the result first found by C. Ewerz [16]. 




a) 



b) c) d) 

Figure 14: Examples of diagrams defining B^$ and B^q. 



Next the //, dd, and 55 terms. We find, in addition to a set of pieces which have the same structure as 
in the photon case, a new structure, Qa-o- In detail: 

B 4;0 (l,2,3,4) = D§ 2} (1,2,3,4) + D^ } (1,2, 3, 4) + £>gj 3} (l, 2,3,4) + Q 4;0 (l, 2, 3,4). (63) 



Here the first three terms are dipole-like, and they follow the reggeization pattern found for the photon 
scattering: 



D^f (1,2,3,4) 



with the color tensor 



_|_ .7 01020403 



D 2 { ^' } (123,4) + D^ } (234, 1) - D^ } (14,23) 



(64) 



Z^ } (124,3)+ Z^ } (134,2) -D^f (12,34) -D^(13, 24)]}, 



r 0102^0304 j 01026 j 60304 
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2iV c 4 4 



(65) 



In the next section we will study the rapidity evolution of these terms, and we will confirm that they follow 
the photon impact factor to all orders. 

The new structure which has no analogue in the case of the photon looks as follows: 



Q4;o(l,2,3,4) 
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[£ 3;0 (12,3,4) +£ 3; o(34,l,2)] + 
[£ 3;0 (13,2,4) + £ 3;0 (24,1,3)] + 
[£ 3; o(14,2,3) + £ 3;0 (23,1,4)]. (66) 
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The function E is the same as in the odderon case, and, in particular, all three quarks participate in the 
interaction. The t-channel state which couples to Qa-$ is Bose symmetric 

Q 4]0 (a(l),a(2),a(3),a(4)) = Q 4 . (l, 2,3,4) (67) 

for any permutation a, and it is gauge invariant: 

Q4-o(k\, ^2,^3,^4) = for any kj — ► 0. (68) 

This property may be proven using the identity for color tensors valid for N c = 3 : 

jai<22b ( jba-ia4 , jaia^b jba-ia^ , ja\a^b jba^az __ _ / r 0102 £0304 1 £ 01 03 r 02 04 , 10104^8203 -i /gg^i 

o 

The analysis in the following section will show that this novel piece of the baryon impact factor couples 
a three-gluon t-channel configuration in which one of the reggeized gluons is an even-signature d-Reggeon. 
The overall C parity therefore is positive. 



5 Integral evolution equations 

In this section we study higher order corrections in the (generalized) leading logarithmic (log s) approxima- 
tion. The all-order sum of these terms will be represented by integral equations [14, 15], written for Mellin 
moments of the multiple discontinuities with respect to the energy s. In our notation the dependence of the 
amplitudes Bn and Bjy (and also Dn, En and Qn) on the Mellin variable u is implicit. 

Let us begin with the C-odd configurations. In the case of three gluons, the impact factor E 3 -q is simply 
replaced by the Green's function function E 3 , which satisfies the BKP equation for three odd signature 
Reggeons, with the initial condition given by E 3 -q: 

L-Y,P(ki)) E 3 = E 3 . + £ K 2 ^ 2 (r,s) ® E 3 , (70) 

V i J (r,s) 

where K 2 ^2 is the real emission part of the BFKL kernel, and the odderon state with the full color structure 
reads 

S3(l,2,3) = d a ^ a *E 3 (l,2,3). (71) 

The four gluon case has been studied in [16], and we simply quote the solution: 

B 4 (l,2,3,4) = (72) 



f aia2b d basa4 E 3 (12,3,4) + f aiasb d ba2a4 E 3 (U,2,4) + f aiaib d ba2Cti E 3 (U,2,3) 
+ f a2a:ib d baiai E 3 {23, 1, 4) + f a ^ b d baiaz E 3 {2A, 1, 3) + f a:iClib d baia2 £ 3 (34, 1, 2) 



where the function E 3 has been defined before in (70). Clearly, the solution is saturated by a reggeizing 
contribution: in each term, one of the three /-Reggeons splits into two elementary gluons. 
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We now turn to the C-even contributions. The integral equations for the multiple discontinuities read 
(up to four gluons): 



to 



U) 



X>( fc * 



£/?(*< 



Y,m 



B 2 



e, 



B A 




(73) 



(74) 



(75) 



The integral kernels driving 2 — > 2,3,4, ... Reggeon transitions, that appear in the above equations include 
the color structure, and they are defined in Ref. [14, 15]. The gluon Regge trajectory (3(k) will be specified 
below. The case of two gluons is the simplest one: Bi satisfies the BFKL equation. According to the 



structure of the inhomogeneous term in Eq. (49), B-i can be written as the sum of three terms D 
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{2,3}, 



1,2)+D^(1,2) + D™{1,2) 
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,{m} 



,{i,j} 



£/3(*i) I>} w = DJ$> + K 2 -, 
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{i,j} 



(76) 



(77) 



i=\ 



In the case of three gluons, B% is given by the sum of three dipole-like components (cf. (60)): 
B 3 (l,2,3) = D< 1>2} (1,2,3) + Dp 1 (1,2, 3) + Df >3} (1,2,3), 



(78) 



where each term consists of three reggeizing pieces: 



d 3 {jj} (i,2,3) = l - g r a2az 



D 



{id}, 



•>{i,j}, 



){»J'}/ 



(79) 



(12, 3) - B\ 4, "(13, 2) + Z^'' J '(23, 1) 

This structure coincides with the photon case. 

The case of B4 is more complex. Following our result for the baryon impact factor in Eq. (63) we 
decompose B4 in the following way: 



,{1,2}/ 



,{1,3} 



{2,3}, 



£4(1,2,3,4) = £>r Z| (l,2,3,4) + 1^(1,2,3,4) + D\^ (1,2,3,4) + Q 4 (l,2,3,4) 
,{m'} 



(80) 



For the dipole-like pieces D\ we make use of the 'reduction procedure' developed for the photon case. 
Namely we decompose each D\ into a reggeizing and an irreducible contributions 

D£ J} (1,2,3,4) = Z^ j};fi (l,2,3,4) + Z)| ij};/ (1,2,3,4), (81) 
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with the reggeizing contribution given by 



D 



{i,j};R t 



1,2,3,4) 



_„2 I Aa\aiaza/± 



,{i,j}. 



,{i,j} 



{i,j}. 



Dl Jt (123, 4) + L>* w (234, 1) - D\' n {U, 23) 



(82) 



+ d 



<2ia2<24<23 



D 2 {ij} (124,3) +^ J} (134,2) - D^' j} (12,34) - D^ j} (13, 24) ] } 



The reggeizing contributions are simple BFKL ladders with one reggeizing gluon splitting into three gluons 
or both reggeized gluons each splitting into two gluons. The irreducible contribution, containing the 2^4 
Reggeon transition vertex, is illustrated in Fig. 15. 




Figure 15: The component D\ 1 ' (1,2,3,4). 

These results provide further evidence that the 'dipole-like' pieces of the baryonic impact factor really 
behave in exactly the same way as the color dipole photon impact factor. In particular, if we would apply 
the large N c limit to the gluon evolution below the impact factor (which, of course, would be inconsistent 
with our finite- N c baryon), the four gluon system below the 2^4 transition vertex would split into two 
non-interacting BFKL ladders, and we would arrive at the first iteration of the BK equation. 

After subtracting, from i?4(l,2,3,4) in (80), these dipole-like contributions of the baryon we are left 
with Qa,. As Qa appears at level of four gluons, its evolution equation has simply the BKP form: 



u> 



£>(** 



a 




(83) 



Making use of the experience with D4, we decompose the amplitude Q4 into a reggeizing piece Q^ and an 
irreducible contribution Q\: 

Q 4 (l,2,3,4) = Qf(l,2,3,4) + $(1,2,3,4). (84) 

Going through steps similar to the ones outlined in [14, 15] we find that the reggeizing pieces Qf preserve 
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Figure 16: The new 3^4 transition vertex W . 



structure of Q^q: 
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[#3(13,2,4) +# 3 (24,1,3)] + 
[# 3 (14,2,3) +# 3 (23,1,4)]. 



(85) 



As seen from the color and momentum structure, the three gluon state coupling to Qa_-$ consists of three 
reggeized gluons, one of which is in a d state and decays into two elementary gluons (the pieces proportional 
to color tensors 55 play the role of subtractions; in particular, they are needed in order to satisfy the Ward 
identities). This state, consisting of two odd signature /-Reggeon and one even signature <i- Reggeon, belongs 
to even C, i.e. to the Pomeron channel. 

The remaining piece, Q A , contains a new transition vertex. We illustrate this contribution in Fig. 16. 
This vertex describes the transition from the three Reggeon state consisting of two / and one d Reggeon to 
four / Reggeons. In more detail, the vertex may be completely decomposed into non-connected pieces of 
two types: (i) the incoming d Reggeon together with one of the / Reggeons makes a transition into three 
/ Reggeons whereas the remaining / Reggeon acts as a (^-channel) spectator, and (ii): two / Reggeons 
interact via the BFKL kernel and the d Reggeon splits into two / Reggeons. The explicit functional form 
of the vertex W, acting on the three Reggeon state fo is the following: 



(W^ 3 )(l,2,3,4) 



-g 



12 3 4 1324 1423 

5 aia2 5 a3a4 {W(f>3) + 5 aias 5 a2a4 (Wcb3) + 5 aia4 5 a2a3 (W(/)3) 



(86) 



where 



1234 

(W<fe) 
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(12)o3 ..4 11o2-(34) 

Q 4>3 + 77 Q 4>3 



+ [3^4] + [l<-»3, 2<->4] + [l<->4, 2<-»3]. 



(87) 



Let us stress that this vertex acts on a completely symmetric function ^3 with three arguments, which 

inherits its structure from #3. This vertex is closely related to a 3 — > 4 vertex found in Ref. [31] in an 

123 
analysis of jet production amplitudes at small x. The symbol Q denotes the integral operator G(l,2,3), 
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introduced first in [14] and further investigated in [32]. It acts on a two gluon amplitude, (j) 2 , and describes 
a transition to three gluons. It consists of two pieces: 

G(l,2,3) = Gi(l,2,3) + G 2 (l,2,3), (88) 

where the first one contains s-channel gluons ('connected part'), the second one reggeizing pieces ('discon- 
nected part'): 

(G 1 (j)2)(kl,k2,k 3 ) = 

d 2 l ((k 2 +k 3 ) 2 l 2 (k 1+ k 2 ) 2 (k-l) 2 k 2 (k-l) 2 l 2 h 2\, (lh n (kQ , 

+ ,u i T^2 77 17YT7T. 1 17vi ~ k I Mh* ~ I), (89) 



and 



(2vr) 3 V (l-h) 2 (k-l-k 3 ) 2 (I - k ± ) 2 (k - I - k 3 ) 

N c g 2 (G 2 <j) 2 )(k 1 ,k 2 ,k 



J ^l 2 (k-l) 2 {[f3(k 2 + k 3 )-f3(k 2 )](27T) 3 5^(l-k 1 ) 

+ [f3(k 1+ k 2 ) - (3(k 2 )}(2irf5W(l-k 3 )} Ml,k-l), (90) 

with the gluon trajectory function 

f d 2 l k 2 1 

P(ki) = -N c g 2 J {27T)3l2 + {k t t _ l)2{kt _ l)2 , (91) 

and k = k\ +k 2 + k 3 . In (87) we have used a short-hand notation for the argument structure introduced in 

123 -4 

Ref. [31]: in the first term, Q <p 3 , <f> 3 is the three gluon amplitude above the vertex W where the rightmost 
Reggeon (momentum k^) is a spectator, and the G operator acts on the two left Reggeons, turning them 

(12)o3 -4 

into the three gluons with momenta ki, k 2 , and k 3 . In the fourth term, Q <p 3 , Reggeon 4 is, again, a 

spectator, and the G operator (with zero momentum in the second outgoing gluon) equals the BFKL kernel 

acting on the two leftmost gluons inside 4> 3 : after this BFKL interaction the leftmost gluon splits into two 

lo2-(34) 

gluons with momenta k± and k 2 , and the other one carries momentum k 3 . Finally, in the last term, Q <p 3 , 
the rightmost spectator now splits into two gluons with momenta k 3 and £4, and the G operator, like in the 
previous term, equals the BFKL operator with outgoing momenta ki and k 2 . 

The full vertex W in (86) is gauge invariant, infra-red finite and Bose symmetric. As the vertex is 
expressed in terms of the function Q, it is also Mobius invariant [32]. Finally, there is no violation of 
signature conservation: the incoming three Reggeon state, consisting of one ci-Reggeon and two /-Reggeons, 
has even signature; the same holds for the outgoing four Reggeon state (four /-Reggeons). 

As a result, the baryonic impact factor introduces a new contribution to the Pomeron channel which has 
no analogue in the photon dipole factor. 
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6 Baryon wave functions in the coordinate space 



The baryon wave function in transverse position space may be easily obtained by the Fourier transform: 
e> A2)Aa ({<*}, {r,},P) = f ^1^^ ¥<**>*» ({«,}, (ft},P) exp fi f ft . r< l • (92) 



2vr 2vr 2vr 



({aJ,{Pi},-P) exp i J^ft -rj 
V i=l / 



The result takes a rather simple form: 
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3 ({a i },{r i },P) = J\f aia 2 «3 exp 



M 2 



E a i ( r i ~ R ) 



exp(iP • R) x 



(93) 



A M 8 X , Ax 5a, -a 2 5 a , a 3 [ (r 2 - il) • ??a ] [ (n - r 3 ) • t?_a ] 



+ A M <5 A , -Ai 5a, a 2 ^a, a 3 [ ( r i ~ R) ■ Vx ] [ ( r 2 - r 3 ) -»|-a] + 



2i 5 a ,a 1 5a,-a 2 5a,-a 3 [(r 2 -r 3 ) ■ rf X ] + 



2«5 a ,-a 1 5a,a 2 5a,-a 3 [(*"i - r 3 ) ■ »7a 



where P denotes the light-cone center of mass position vector, 



OLiTi. 



i=\ 



(94) 



The form of the wave function given by Eq. (93) which follows from the Ioffe current shows in detail the 
angular momentum structure of the baryon and the correlations between the angular momenta and quark 
helicities. In particular, each scalar product of the type (fi — R) r) X clearly indicates a rotation of quark 1 
around the baryon center-of-mass with the orbital angular momentum z-component, L z , equal to A. Terms 
of the type {t\ — J" 3 ) r) X correspond to a similar rotation within the quark pair (1,3). Thus, in the massless 
quark case, all components of the baryon carry a non-zero angular momentum L z for the Ioffe operator. An 
inspection of the momentum space expressions (32) shows that for the massive quark case, one may have 
Ioffe baryon wave function components with L z = 0. 

Using Eqs. (43) and (39), one may express the baryon impact factors B^ ({^};P,P') via the overlap 
function T AA ({li};P,P r ) defined in the coordinate space: 



E j{d 2 ri][d> 

Ai,A 2 ,A 3 
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(95) 



exp -i^lj'fj 
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a ({«,}, {r,};P). 



It follows from Eqs. (45), (46) and (95) that the normalization condition for the wave function reads: 
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({a,}, {r,};P) = 5^{P-P'). (96) 
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It is instructive to evaluate a contribution to the baryon two-gluon impact factor [<5£>2;o] corresponding 

to a dipole-like piece, e.g. to D^.q , in the coordinate representation. The gluon color labels are a\ and a 2 
and momenta are denoted by k\ and k 2 respectively. One obtains: 



5B^({h};P,P') 



,(Ai,A 2 )A 3 



{1,2} 1 S aia2 
= - i-iq) x 

2 v yj 2N r 



(97) 



x Y, l[d 2 ri}[dai] tf^>* ({<*}, -fa}; J") 
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- ifci-ri — ifc2-i , 2 



- jfci-r2 — iA^-ri 



Assuming, for simplicity, the forward kinematics, k\ = k = —k 2 , one may rewrite the eikonal factors in the 
last line of (97) in a factorized form, found in the case of the color dipole scattering, 



1 



Jk-(r 2 -r 1 ) 



1 



= ifc-(r 2 -ri) 



(98) 



This equivalence of the structures holds also beyond the forward limit (note that, for nonzero P, P' the 
wave functions & \}' 3 contain the phase factors e ). In Eq. (97), the prefactor 1/2 in the first line 
reflects the relative weight between the color dipole scattering amplitude and the scattering amplitude of 
the dipole-like components of the baryon. 

7 The quark— diquark limit 



In many phenomenological applications the nucleon is represented as a bound state of quark and a tightly 
bound diquark. The transverse size of the diquark is then assumed to be much smaller than the size of 
the baryon, and the diquark state emerges in an anti-triplet color representation. In this approximation 
the baryon should resemble an (asymmetric) color dipole. It is interesting to analyze the properties of 
our baryon impact factor in this limit. Formally, the quark-diquark limit corresponds to the limit where 
the transverse separation of two quark lines shrinks to zero, and a t-channel gluon no longer distinguishes 
between the two quark lines. In momentum space, as seen in (95), the overlap function then only depends 
upon the sum of the momenta of all gluons coupled to the two coinciding quark lines. To be definite, let us 
assume that quarks 2 and 3 move close to each other. Then all overlap functions F degenerate to a function 
^? 1 (23) w jth on ly t wo arguments: 



F(k 1 ,k 2 ,k 3 



3^2 



F^ikukz+k 



3 ■ 



(99) 



(note that the limit F^ 2 ^(ki, k 2 ) is not necessarily symmetric in its arguments). Applying this argument to 
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the three dipole-like terms in (49) we immediately see that the dipole-like component D\.^ vanishes if lines 
2 and 3 are contracted: this is the well-known limit of a dipole with vanishing size (color transparency). In 
more detail, (52) shows that all terms in this impact factor tend to i ?1 '- 23 '(0,fci + k 2 ), and they cancel due 



to opposite signs. The remaining dipole-like components D^-rj and D^.q become equal: 
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with 

D^ 23)} (k 1 ,k 2 ) = ^ [F 1 ( 23 )(0,fc 1 +fc 2 ) + F 1 ( 23 )(fc 1 +A ;2 ,0)-F 1 ( 23 )(fc 1 ,fc 2 )-F 1 ( 23 )(A: 2 ,fc 1 )]. (101) 

As we already discussed at the end of Section 4.2, in (49) each dipole-like term carries a color factor 1/2, 
compared to a genuine color dipole factor seen in a color singlet quark-antiquark system. Since in the 

12 3 \ 11 2\ {1 3^ 

quark-diquark limit -D 2 .q vanishes and the contributions from D^.q and D^.q coincide, this part of the 

Ji (0 3 )\ 

baryonic impact factor adds up to a standard dipole contribution D 2; o(ki,k 2 ) = 2Z?2-o (ki,k 2 ). 

Next, we turn to the three-gluon impact factors. In the Pomeron channel, one finds only reggeizing pieces 
of the quark-diquark dipole impact factor. In the odderon channel, the function E(ki,k 2 , k 3 ) degenerates to 
the structure found in the 7* — ► r/ c transition impact factor, which couples only to the Bartels-Lipatov-Vacca 
(BLV) odderon [33] but not to the Janik-Wosiek solution [34]: 

E 3 .,o(k 1 ,k 2 ,k 3 ) ^ E^ 23)} (k 1 ,k 2 ,k 3 ), (102) 

with 

E^WfoMM) = ^ |> (23) (*i,*2+*3) - F 1 ^(k 2 +k 3 ,k 1 ) + 

+ F 1 ^ 23 \k 2 ,k 1 + k 3 ) - F 1 ( 23 )(fci+/b 3 ,fc 2 ) + F 1 W(k 3 ,k 1 +k 2 ) - F 1 W(k 1 +k 2 ,k 3 ) + 

+ F 1 W(k 1 +k 2 + k 3 ,0) - F 1 ^(0,k 1 +k2 + k 3 )\ . (103) 

For the four gluon case, one finds the standard reggeizing pattern of -D 2 .q and of E^.q in the 

Pomeron and the odderon channel, respectively. The structure Q^-q vanishes in the quark-diquark limit. 
This is the result of a nontrivial cancellation of all three lines of Eq. (66), making use of the identity (69). 
The pattern given by the impact factors in the small diquark limit is preserved by the small x evolution, in 
particular Q4 vanishes. 

In summary, we have verified that, in the quark-diquark limit, the baryon reduces to a dipole-like 
object with an asymmetric wave function, as it was expected. Conversely, our analysis shows that, outside 
the diquark limit, the baryon impact factor contains a new piece (related to Q^-o) which accompanies the 

{2 3| 

appearance of the third dipole-like term, D 2 .q . A more detailed study of the question, to what extent the 
baryon wave functions actually favors a diquark state, should start from the Fourier transform of the overlap 
function, (95), which describes the distribution of the quarks in transverse coordinate space. Further work 
along these lines is in progress. 



8 Discussion 

In this paper we have investigated the high energy behavior of a baryonic state. We have studied the 
structure of a baryonic impact factor, its coupling to multi-gluon exchanges and the rapidity evolution of 
the ^-channel gluon states. We found it convenient to follow very much the same approach, which has been 
developed and used for the high energy behavior of a virtual photon (or a heavy quarkonium state). For 
the scattering of such mesonic states, in the leading logarithmic approximation and in the large-iV c limit, 
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the high energy behavior allows for the interpretation in terms of color dipoles, and one of the motivations 
of our investigation was the question to what extent this attractive physical picture can be used also for the 
scattering of baryonic states. 

Compared to the quark-antiquark system created by the photon (or a heavy vector meson), the high en- 
ergy scattering of baryonic systems consisting of three quarks shows similarities, but also striking differences. 
First, there is a component of the baryonic impact factor in which two of the three quarks interact with the 
target whereas the third one acts as a spectator. Here the two-quark subsystem behaves very much in the 
same way as the color singlet dipole of the quark-antiquark system. In particular, the rapidity evolution 
is the same as in the case of a virtual photon. This configuration, however, extends beyond the picture of 
a small "diquark state": we have shown that, in the diquark limit, we recover the dipole picture. But the 
spectator quark is not necessarily linked (in transverse space) to one of the participating quarks, and our 
analysis includes also this more general configuration. Second, there is the piece of the baryon impact factor 
to which the C-odd three gluon state (odderon) couples. Third, a new piece of the baryonic impact factor 
exists which couples to a C-even three gluon t-channel state, and there is a new vertex which describes the 
transition from this three gluon state to the four gluon (two Pomeron) state. In the quark-antiquark case, 
there is no analogue of this contribution. 

This third piece may actually be quite essential for the restauration of s-channel unitarity in baryon 
scattering and can therefore not be neglected. Namely, let us consider the scattering of a hypothetical heavy 
baryon on a large nuclear target; this represents the analogue of the Balitsky-Kovchegov problem for the 
color dipole scattering. Based on our results, the baryon scattering amplitude B can be written symbolically 
as a sum of the following pieces: 

C— even 



,{1,2} . „{1,3} . ^{2,3} 



C*-odd 



B = Vp^ + Vf'^ + V\^ + Q A + £ 3 ■ (104) 

Here the first three terms, T)\ , stand for the dipole-like contributions in which the baryon couples to the 
same two-point gluon correlator as the color dipole in the scattering of a virtual photon. The strength of 
this coupling, however, is only 1/2 of that for the photon dipole. The pieces Q4 and £3 probe three-point 
gluon correlators: the C-even and C-odd ones respectively. As it was observed in the case of the color dipole 
in deep inelastic scattering, where only a single BFKL Pomeron could couple to the dipole, we again see no 
indications of a direct two Pomeron coupling to the valence degrees of freedom of the baryon. If we assume 
that the two-gluon distribution probed by the first three terms in Eq. (104) is consistent with saturation 
of the black disc limit for color dipoles of the sizes given by the baryon geometry, then the T-matrices for 
each of the Pg components would tend to 1/2, and the total contribution of the dipole-like pieces to the 
baryon T-matrix would amount to 3/2. This would mean that s-channel unitarity can be maintained only 
if Qa and £3 give a combined contribution to the T-matrix smaller than —1/2. Thus, the three- Reggeon 
states Q4 and £3 seem to be essential to guarantee the s-channel unitarity. Interestingly enough, one 
might go even further and arrive at a quantitative prediction: if one postulates that the T-matrices, both 
for the color dipole and the baryon scattering at very large energies saturate the unitarity limit — one then 
finds that in the black disc limit: (i) the C-odd three point gluon correlator should vanish; this comes from 
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the requirement that both proton and anti-proton scattering should reach the black disc limit, despite the 
fact that the amplitude £3 has opposite signs in these two cases; (ii) the C-even three point correlator is 
strongly constrained: when coupled to the impact factor Q^-o it must lead to the scattering amplitude equal 
to —1/2. In the diquark limit, both D\ ' and Q4 vanish, and unitarization proceeds in the same way as in 
the dipole case. 

We interpret these results as a strong indication that, in the context of baryon scattering, QCD Reggeon 
field theory has to be extended beyond the theory of BFKL Pomerons and their interactions. First, it is 
difficult to justify the large- N c limit, which, in the scattering of virtual photon and mesonic states, allows to 
reduce the evolution of BKP states consisting of 2n-gluon to the propagation of n BFKL Pomerons. Second, 
the three gluon state (and its BKP evolution) seems to play an important role, not only in the odderon 
channel. As we have pointed out, this phenomenon is closely connected with the existence of the d Reggeon, 
the even signature partner of the (odd signature) reggeized gluon. 

On a deeper level one may speculate that there exists an intimate connection between the number of 
valence objects in the impact factor in the fundamental SU(N C ) representation and the maximal number of 
Reggeons in the BKP state which couple to the impact factor. For the quark-antiquark color dipole only 
the two-Reggeon BFKL Pomeron couples, and for the baryon containing three quarks we have both two- 
and three-Reggeon states. We may conjecture that the number of the different BKP states that couple 
to the baryon in SU(N C ) gauge theory is related to the number of Casimir operators of the gauge group. 
There exist two Casimir operators of the SU(3) gauge group, and QCD Reggeon field theory (whose basic 
degrees of freedom are the reggeized gluons) exhibits two 'fundamental excitations' which, in the leading- log 
approximation, are represented by the two-gluon BFKL Pomeron and by the three-gluon odderon state. For 
a high energy SU(N C ) baryon we expect that the impact factor, consisting of N c quarks in the fundamental 
representation, would exhibit all the 2, 3, . . . , N c gluon states, and it would hint that the number of 
fundamental glue excitations may be related to the N c — 1 Casimir operators of SU(N C ). It seems natural 
that the gauge group invariants should be mapped onto gauge invariant BKP states. The explicit connection, 
however, has not been yet established. 

Turning to more practical and phenomenological applications, in this paper we have considered a baryonic 
state consisting of three massive quarks being in a proton-like configuration. One can view such a 'heavy 
baryonium' state as a convenient theoretical laboratory, very much in the same spirit as previous work on 
high energy QCD has made use of 'heavy onium' states. On the other hand, we feel that our results might 
also allow for immediate phenomenological applications. In particular, we have proposed a relativistic 
invariant model of the proton wave function, including the helicity structure and correlations between 
helicities and quark angular momenta. Both the model itself and the calculational technique applied may 
be useful in studies of polarized scattering of the proton and of the proton form-factors. Another potential 
place of interest is the intermediate t region of proton-proton elastic scattering where, in the days of ISR 
experiments, a very simple three gluon model had a striking phenomenological success [35]. It should also 
be quite interesting to study other applications of the model in the context of elastic pp and pp scattering 
and exclusive diffraction at RHIC, Tevatron and the LHC. Finally, we would like to view our study as a 
preparation for a QCD analysis of multiple scattering in pp collision at the LHC. 
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A Appendix 

A.l Spinorial matrix elements 

The calculations of the baryon wave functions and of the baryon scattering amplitudes are performed using 
the light-cone formalism summarized in [23]. 

Thus we employ the spinor basis defined by 




1 



and 



where 
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(105) 
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(107) 



in the Dirac representation, and the Dirac matrices J3 and a are related to the 7-matrices through J3 = 7 
and a s = 7°7 S ; m is the mass of a fermion (or an anti-fermion) . In the infinite momentum frame, when 
p + — ► 00 these spinors tend to the helicity eigenstates, u-\i(p) — ► u±(p), v-n(p) — ► v±(p). 

In the calculation of the baryon — > quarks transition amplitudes it is sufficient to employ spinor matrix 
elements given in the following tables. Note that we consider a general case in which the masses of the 
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spinors u (or v) and u' are given by m and m', respectively. 



Matrix clement 
u' y {p)...u x {q) 



A -»V 

T -» T 
I -» I 



A -»V 

T -» I 
I -> T 



"'(p) 7 + "(g) 



"'(p) -,- "(g) 
a/p+ ' v^ 



^M^s "(g) 



^4t[(p-»?t)(9-^±) + mm 7 ] 



r&^F +*£*& 



P T <3 n 
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±V± 



aim m 



Matrix element 
v' x ,(p)...u x (q) 



A ->A' 

T -» T 
I -» I 



A ^A' 

T -* I 
I -» T 



^'(p) „,+ »(g) 

/„-l_ ' / L 
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Vp+ ' a/?"*" 
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^qr[(p-r?q=)(9-»7±) ~ "im'] 



„s P'Ht 1 s q-*?± 

»7± p+ + ^T 0+ 



As an example, we apply the above formulae to evaluate 

[d\ [i {P3)7fiW\(P)] ■ [u\ 1 {pi)j' I Vx 2 (p2)} _ 1 [d\{P3)j + W X {P)] ■ [u Xl (pi)-f~Vx 2 {p2)_ 



p+ pf p+ p+ 



p+ p+ p+ p+ 

1 [d\{P3)j~W \{P)] ■ [^A!(Pl)7 + ^A 2 (P2)" 



'P+p+p+p+ 

[ d\ (ps ) 7± ^a ( P)] • [u\ 1 (pi)Y±vx 2 (P2) 

P+ p+ p\ p\ 



(108) 



for A = Ai = — A2 = A3 = +1. The prefactors: 1/2, 1/2 and —1 on the r.h.s. are the only non- vanishing 
elements of the covariant metric tensor g„ v in the light-cone coordinates. In the calculations we find it 
useful to make use of the following identities for transverse complex vectors r)±: rf + = rf~ , r)\ = 0, 
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r)± rj^z = \r)±\ 2 = 2. Thus, assuming that the light quark masses vanish, we obtain: 



[ ^A 3 (P3 ) 7m ^A ( P)] ■ [v\ 2 {P2)Yux 1 {pi, 

p+piptpt 



2[(P2-*7-)(pi-i7+)]* , 2(p 3 • V-)(P ■ »7h 
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P2P3 

p 2 • *?+ P • V-\ 



Pi 



P^ 



P+P+ 

Pi V- P3 ■ *? 



Pi 



P,3 



(109) 



Using an identity 3 dA 3 (P3)75 = ^3^A 3 (P3), and relation (10), one obtains one of the matrix elements 
described by (21). The matrix elements for all remaining choices of helicities can be derived in the same 
way. 

A. 2 A reduction formula for spinors in high energy limit 

We shall prove the following identity for massive Dirac spinors: 

u(p) q (p + m + k) = 2p ■ q u(p + k) + . . . , (HO) 

which holds, at the leading accuracy in s ~ 2p ■ q, in the high energy limit: s 3> q 2 , k 2 , m 2 , p • k, q • k etc., 
and for k± S> k + ,k~ . This identity is a useful tool for deriving quark scattering amplitudes by multi-gluon 
couplings in the eikonal approximation. Using the spinor equation of motion, u(p)(p — m) = , we get 



u(p) q(p + m + k) — u{p) [ 2p • q + q k ) ~ s u(p) ( 1 + — 



q, k 



fill) 



where we used the fact that the anticommutator {q , k} = 2k ■ q <C s. Furthermore, using the light-cone 
variables, as defined in Sec. 3, we have 



A- 



-2ia a/ 3q a k^ ~ —2ia^ r q kj_, 



where r is the Lorentz index of the transverse coordinates. Thus one obtains 

2 (J rp rZ I 



u{p) q{p + m + k) ~ s u{p) I 1 



pt 



(112) 



(113) 



The matrices a a ' are proportional to the generators of the Lorentz transformations of the Dirac spinors: 



, a l 3 \ „(n\ - -,,(\(,.A^ iifrA o™ ( -rr ^,., a > 3 



exp ( --<J a pu ali \ u{p) = u(A(w)p), u{p) exp I -a a puj ali ) = u(A(w)p), (114) 



where 



(A(io)pr = [A(w)] V, A(w) = exp Qu^O , 
and the generators of Lorentz transformations in the vector representation read 

[^Fj, = 9^ a 9uf3 ~ 9^/3 Qua- 



(115) 



(116) 



3 For a non-zero quark mass m, the relation holds approximately in the large energy limit, d\ 3 (p 3 )y 5 = A 3 gJa 3 (p 3 ) + 

0(m/p+). 
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Since the parameter multiplying <r_ r in Eq. (113) is small, f3 r = k^/p^ <C 1, one may write 

u{p) (l-ia- r (3 r ) = u{p) [expi-i&orF/V) exp(ia 3 r/3 r /2)] + 0(f3 2 ), (117) 

where we used the identity 7_ = 5(70 — 73)- This equation corresponds to two subsequent infinitesimal 
Lorentz transformations acting on u{p) with the parameters lu\° = — u^ r = (3 r and uj^' = —^'2 = P r 
(and all other components ui^ = 0)- This is an infinitesimal boost along the transverse direction /?, and an 
infinitesimal rotation in the plane spanned by the transverse vector around the z-axis. Using (114) one sees 
that, in leading order in f3 r , the boost transforms p in the following way : p° — ► p°, p — ► p + p°0, p 3 —> p 3 , 
and the rotation acts as: p° — > p°, p — ► p + p 3 /?, p 3 — ► p 3 - Thus one obtains 

u{p) (l-ia- r (3 r ) = u(p') + 0(/3 2 ), (118) 

with p' = (p°,p + /9p + ,p 3 ). This proves Eq. (110). The equation for multiple eikonal couplings, Eq. (15), 
follows immediately from Eq. (110), after all spinor contractions, q (p — h\ — . . . — ki + m) q — 2p ■ q q, are 
executed. 
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